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Field of the Invention 

The present invention relates generally to the field of digital signal processing, 
10 and more particularly to the use of Gaussian pulses to perform Time Domain 
Reflectometry to locate and characterize faults in a Device Under Test. 

Description of the Related Art 

In the field of digital signal processing there is often a need to detect a signal pulse 
15 in a transmission medium, such as an electrical cable or optical fiber, or an electronic or 
optical device. For example, in electronic or optical systems, when a signal encounters a 
change in impedance in the transmission medium, a partial reflectance of the signal 
occurs, decreasing the amplitude of the propagated signal. These impedance changes can 
substantially degrade the performance of a device or transmission medium, and so a great 
20 deal of effort has been spent detecting and analyzing such interconnect discontinuities. 

One technique which has been used to great effect in this area is Time Domain 
Reflectometry (TDR). TDR involves measuring reflection in an unknown device or 
medium in comparison to reflection due to a standard impedance. In other words, TDR 
compares reflected energy to incident energy on a single-line transmission system. 
25 Known incident stimulus is applied to the standard impedance and intentionally 
propagated toward the unknown device or medium. Reflections from the unknown 
device or medium then propagate back to the source, and amplitude and time of the 
reflected signal(s) is compared to the incident stimulus. The reflected signal magnitude 
(and possibly the waveform) is a function of the incident signal magnitude and the nature 
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of the impedance change. The time elapsed between the detection of the incident and the 
reflected signal is a function of the overall distance traveled and the velocity of 
propagation of the signal. Thus, by detecting and timing the reflected pulse with respect 
to the original pulse, the distance to an interconnect discontinuity may be determined. In 
5 this way, for example, a flaw, such as a kink, may be located precisely in the device or 
medium. 

Fast sampling systems with TDR capability can be very useful for studying 
interconnect discontinuities. In a typical TDR system, a step generator, such as a 
switched current source, is used to generate a step pulse, which is propagated through a 
10 T-connector to a digitizer and the device under test (DUT) respectively. The reflected 
signal (from the DUT) is received by the digitizer and the superpositioned waveforms 
analyzed to determine the elapsed time, and to compare amplitude and waveform of the 
pulses. 

However, there are a number of drawbacks to this approach. The use of a step 
15 function or signal is problematic because the steep rise of the leading edge of the step 
entails many high frequency signal components. In other words, the step signal is a wide 
bandwidth pulse. Transmission media and electronic and optical components behave 
differently for different frequencies, leading to such distorting effects as dispersion and 
dissipation of the pulse energy. The effect of dispersion is particularly troublesome for 
20 step pulses in TDR applications because the leading edge of the step is used to determine 
the timing of the pulse. The dispersion of high frequency components in the edge can 
smooth the edge to such a degree that precise timing of the reflected pulse may be 
difficult or even impossible. 

Current systems developed to perform TDR with step pulses involve extremely 
25 high speed sampling and tend to be very expensive, costing on the order of one hundred 
thousand dollars each. The use of modulated Gaussian pulses could avoid many of the 
problems associated with step pulses in that Gaussian pulses are extremely smooth, have 
low bandwidth, and are particularly suitable for closed form analytic operations and 
representations. However, there are currently no known systems or methods to perform 
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TDR with Gaussian pulses, due to difficulties in accurately detecting, timing, and 
characterizing Gaussian pulses in a noisy medium. 

Therefore, systems and methods are desired to detect and characterize modulated 
Gaussian pulses in a noisy medium. 
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Summary of the Invention 



A system and method for detecting and characterizing Gaussian pulses is 
presented. A system and method for performing Time Domain Reflectometry (TDR) 
using Gaussian pulses is also described. The system may include a computer, comprising 
a CPU and a memory, wherein the memory is operable to store one or more software 
programs for performing TDR, and wherein the CPU is operable to execute the software 
programs. The system may also include an arbitrary waveform generator (AWG) coupled 
to the computer, and a digitizer coupled to the computer and the AWG. The system may 
be coupled to a Device or Medium, hereafter referred to as a Device Under Test (DUT). 
It should be noted that the DUT may be any kind of device, including a stand alone 
device, a PC board, an instrument, an electric or optical circuit, or an electronic or optical 
transmission medium, among others. 

In the preferred embodiment, the computer system comprises a PCI extensions 
for Instrumentation (PXI) system which includes one or more PXI computer boards or 
cards plugged into a PXI backplane, such as a "PC on a card", housed in a PXI chassis. 
In other words, the PXI cards may comprise the memory and CPU which are operable to 
respectively store and execute one or more computer software programs implementing 
the present invention. The PXI system may also include or couple to a display, such as a 
monitor, for displaying visual information, such as results, to a user, as well as an I/O 
interface for receiving input and sending output to external systems or components. In 
one embodiment, the display may be comprised in the PXI chassis. In another 
embodiment, the display may be external to the PXI chassis. In one embodiment, the I/O 
interface may also be comprised on a PXI card. 

The AWG may be operable to generate a Gaussian pulse and transmit the 
Gaussian pulse to the digitizer and to the Device Under Test (DUT). The DUT may be 
further operable to reflect at least a portion of the transmitted Gaussian pulse in the form 
of one or more reflected pulses to the digitizer. A digital signal containing a modulated 
Gaussian pulse and signal noise may be received, such as from the digitizer, from some 
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other external system, or from the memory medium of the computer system. The signal 
may also include one or more reflected pulses. For example, in the TDR system 
described above, the AWG may generate a Gaussian pulse and transmit the Gaussian 
pulse to the digitizer and the Device Under Test (DUT). The DUT may reflect at least a 
portion of the transmitted Gaussian pulse to the digitizer, comprising the one or more 
reflected Gaussian pulses. The digitizer may receive and digitize the signal comprising 
the transmitted Gaussian pulse and the one or more reflected Gaussian pulses, and store 
the digitized signal in the memory. 

An estimation of N Gaussian pulse parameters may be determined for a Gaussian 
pulse comprised in the signal, where N is greater than or equal to one. The Gaussian 
pulse parameters may include Op (the inverse of the Gaussian pulse variance), t p (the 
Gaussian pulse time shift), and C0p (the Gaussian pulse carrier frequency or base 
frequency). The pulse variance is related to the width of the Gaussian pulse. The pulse 
time shift is a measure of the time interval from some arbitrary origin to the peak of the 
Gaussian pulse. The pulse carrier frequency refers to the frequency of the modulation of 
the Gaussian pulse. In one embodiment, the estimation of Gaussian pulse parameters 
may be determined using any of various prior art techniques, such as the zoom-in 
approach, for example, which is well known in the art. 

Once the estimation of the Gaussian pulse parameters is determined, a plurality of 
permutations of the estimated Gaussian Pulse parameters may be generated by adding or 
subtracting a small delta amount to or from each parameter value, producing a plurality of 
parameter sets representing a corresponding plurality of permuted Gaussian pulses or 
waveforms. If there are N estimated parameters, then M permutations of the N 
parameters may be generated, where M is greater than or equal to N-L The original 
coarse estimated parameters and the M permutations comprise M+l sets of parameters, 
such that for M = N-l, the minimal number of parameter sets (M+l) = N may be used to 
determine refined values for N parameter variables, described below. It should be noted 
that in some embodiments, the coarse estimation may not be used in determining the 
refined values. 
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In the preferred embodiment, M is greater than or equal to N so as to allow 
overdetermination of the N parameters. For example, where the pulse parameters are Op, 
t p , and Op , N = 3, and so M = 3 permutation sets may be generated. It has been 
determined that, due to diminishing returns of increased overdetermination, in the 
preferred embodiment, letting M = N provides for the greatest benefit in refining the N 
estimated parameters. 

In one embodiment, an inner or dot product, Qi, may be calculated between the 
received signal and each of the M permutations of the estimated Gaussian Pulse 
parameters, i.e., between the received signal and each of the M permuted Gaussian pulses 
or waveforms. This calculation is facilitated by the fact that for Gaussian pulses there is a 
closed form, or analytic, solution for the inner product. Said another way, each of the 
plurality of parameter sets may be used in the equation for Q, providing a corresponding 
plurality of inner product equations. The closed form equation of the inner product 
includes terms for the estimation parameters, as well as variables for the as yet 
undetermined parameters of the received Gaussian pulse. 

The plurality of M inner products, Q, may be used to generate a corresponding 
plurality of M linear equations, each of which is a function of a respective one of the M 
estimation permutations and corresponding N parameter variables of the Gaussian pulse. 

The plurality of linear equations may then be solved to overdetermine the 
parameters of the Gaussian pulse. It should be noted that in the case of three parameters 
or variables, a minimum of three linear equations could be solved to determine the 
variable values. However, by overdetermining the variable values through the use of 
more than the minimum number of linear equations, an "averaging" effect results, 
providing more accurate values for the parameters. Using the (over) determined 
parameters of the received Gaussian pulse, various analyses may be performed. For 
example, in the TDR system, a transmitted Gaussian pulse and one or more reflected 
pulses (from the DUT) may each be analyzed to characterize an interconnect discontinuity 
in the DUT. In one embodiment, the parameters may be used to determine the nature and 
magnitude of the discontinuity or the location of the discontinuity. For example, the 
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characterized Gaussian pulses may be used to determine the transfer function of the 
discontinuity, or to calculate the impedance. As another example, the value of the time 
shift parameter t p for each pulse may be used to locate the interconnect discontinuity or 
flaw in the device or transmission medium. 
5 In one embodiment, the method described above may be applied iteratively to a 

signal containing multiple Gaussian pulses, such as may be produced by the TDR system. 
In this embodiment, each time a pulse is detected and characterized, the characterized 
pulse may be subtracted from the signal, leaving a residue containing any remaining 
pulses. This residue may then be used as the input signal, and the next pulse detected and 

10 characterized. Again, the characterized pulse may be subtracted from the signal, leaving 
another residual signal, and so on, until all pulses of interested have been detected and 
characterized. Useful analyses may then be conducted on the resulting characterized 
pulses to extract useful information regarding the DUT, such as the number and nature of 
connection discontinuities therein. 

15 Thus, by using various embodiments of the above-described method, a Gaussian 

pulse in a noisy signal may be detected and characterized. In one embodiment, the 
method may be used to distinguish between an original pulse and a plurality of reflected 
pulses in a TDR system, and to characterize connection discontinuities in the DUT or 
medium. 

20 
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Brief Description of the Drawings 

A better understanding of the present invention can be obtained when the 
following detailed description of the preferred embodiment is considered in conjunction 
with the following drawings, in which: 
5 Figure 1 illustrates one embodiment of a system suitable for performing TDR with 

Gaussian pulses; 

Figure 2 is a flowchart of a method for detecting and characterizing a Gaussian 
pulse, according to one embodiment; 

Figure 3 is a flowchart of a method for determining an estimation of Gaussian pulse 
10 parameters, according to one embodiment; 

Figures 4A-4D illustrate example signals with Gaussian pulses, according to one 
embodiment; 

Figure 5 is a flowchart of a method for performing Time Domain Reflectometry 
using Gaussian pulses, according to one embodiment; 
15 Figures 6A-31B illustrate an example implementation of a system for detecting 

and characterizing Gaussian pulses in a received signal, according to one embodiment; 

Figure 32 illustrates a hierarchy of Virtual Instruments (Vis) included in an 
Adaptive Gabor Expansion VI of the system implementation of Figures 6A-31B; and 

Figure 33 illustrates a hierarchy of Vis included in an Adaptive Gabor Transform 
20 VI of the system implementation of Figures 6A-3 IB. 

While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof are shown by way of example in the drawings and are 
herein described in detail. It should be understood, however, that the drawings and 
25 detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the intention is to cover all modifications, equivalents and 
alternatives falling within the spirit and scope of the present invention as defined by the 
appended claims. 
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Detailed Description of the Embodiments 

Figure 1 - A Gaussian Pulse TDR System 

Figure 1 is a block diagram of a Gaussian pulse TDR system, according to one 

5 embodiment. As Figure 1 shows, a computer system 102 may comprise a CPU 101, a 
memory 103, a display 107, and an I/O interface 105. In the preferred embodiment, the 
computer system 102 comprises a PCI extensions for Instrumentation (PXI) system 
which includes one or more PXI computer boards or cards plugged into a PXI backplane, 
such as a "PC on a card", housed in a PXI chassis. In other words, the PXI cards may 

10 comprise the memory 103 and CPU 101 which are operable to respectively store and 
execute one or more computer software programs implementing the present invention. 
As Figure 1 also shows, in one embodiment, the PXI system 102 may also include display 
107, such as a monitor, for displaying visual information, such as results, to a user, as 
well as I/O interface 105 for receiving input and sending output to external systems or 

15 components. In one embodiment, the display 107 may be comprised in the PXI chassis. 
In another embodiment, the display may be external to the PXI chassis. In one 
embodiment, the I/O interface 105 may be comprised on a PXI card. 

As Figure 1 shows, the PXI cards may also include a waveform generator, e.g., an 
arbitrary waveform generator (AWG) 104, and a digitizer 106 which may each be coupled 

20 to a device under test (DUT) 110. As used herein the term "DUT" includes any of 
various devices that may be tested. The "DUT" also includes any of various types of 
transmission mediums that may be tested, such as electrical or optical transmission 
mediums. Thus, the DUT 110 may be any kind of device, including a stand alone device, 
a PC board, an instrument, an electric or optical circuit, or an electronic or optical 

25 transmission medium, among others. For example, in one embodiment, the AWG 104 
and the digitizer 106 may couple via transmission media 1 1 1 A and 1 1 IB respectively to a 
T-Connector 108, as shown. The T-Connector may be operable to further couple to the 
device under test (DUT) 110 via transmission medium 111C. The AWG 104 may be 
operable to generate a signal pulse which may be transmitted to the digitizer 106 and the 
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DUT 110. The signal pulse may be at least partially reflected by the DUT 110 and both 
the original signal pulse and the reflected pulse received and digitized by the digitizer 106 
for analysis. 

Further details of the analysis and characterization of the signal pulses are 
5 presented below with reference to Figures 2 and 3. 

Figure 2 - Flowchart of a Method for Gaussian Pulse Detection and Characterization 

Figure 2 is a flowchart of one embodiment of a method for detection and 
characterization of a modulated Gaussian pulse. It should be noted that in various 
10 embodiments, various of the steps presented below may occur in a different order than 
shown, or may be omitted, as desired. In other embodiments, additional steps may be 
*f included, as well. 

CO As Figure 2 shows, in 202 a digital signal containing a modulated Gaussian pulse 

m and signal noise may be received. The signal may be received from the digitizer 106 of 

z: 15 Figure 1, or, alternatively, may be received from some other external system, or from the 
If! memory medium 103 of the computer system 102. When the signal is received from the 

O memory medium 103, the signal is presumed to have been previously captured and stored 

r p in the memory medium. 

[3 In one embodiment, the signal may also include one or more reflected pulses. For 

U 20 example, in the TDR system described above with reference to Figure 1, the AWG 104 

may generate a Gaussian pulse and transmit the Gaussian pulse to the digitizer 106 and 
the Device Under Test (DUT) 110. The DUT 110 may reflect at least a portion of the 
transmitted Gaussian pulse to the digitizer 106, comprising the one or more reflected 
Gaussian pulses. The digitizer 106 may receive and digitize the signal comprising the 
25 transmitted Gaussian pulse and the one or more reflected Gaussian pulses, and store the 
digitized signal in the memory 103. 

According to one embodiment, a coarse estimation may be determined for at least 
one Gaussian pulse, as described in 204, then a refined estimate may be determined using 
the coarse estimate, described in 206-212, below. 
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In 204, an estimation of N Gaussian pulse parameters may be determined for at 
least one Gaussian pulse comprised in the signal, where N is greater than or equal to one. 
In the preferred embodiment, the Gaussian pulse parameters may include 0Cp (the inverse 
of the Gaussian pulse variance), t p (the Gaussian pulse time shift), and (Dp (the Gaussian 
5 pulse carrier frequency or base frequency). In one embodiment, the Gaussian pulse 
parameters may further include other parameters, such as the pulse amplitude or phase, 
for example. The pulse variance is related to the width of the Gaussian pulse. The pulse 
time shift is a measure of the time interval from some arbitrary origin to the peak of the 
Gaussian pulse. The pulse carrier frequency refers to the frequency of the modulation of 

10 the Gaussian pulse. In one embodiment, the estimation of Gaussian pulse parameters 
may be determined using the zoom-in approach, which is well known in the art. Further 
details of this approach are presented below with reference to Figure 3. 

Once the estimation of the Gaussian pulse parameters is determined for the 
Gaussian pulse, then in 206 a plurality of permutations of the estimated Gaussian Pulse 

15 parameters may be generated. In other words, for each parameter value, corresponding 
additional values may be generated by adding or subtracting a small delta amount to or 
from the value, producing a plurality of parameter sets representing a corresponding 
plurality of permuted Gaussian pulses. Details of this process are described below in the 
Theory section. In one embodiment, if there are N estimated parameters, then M 

20 permutations of the N parameters may be generated, where M is greater than or equal to 
N-l. In one embodiment, the original coarse estimated parameters and the M 
permutations comprise M+l sets of parameters, such that for M = N-l, the minimal 
number of parameter sets (M+l) = N may be used to determine refined values for N 
parameter variables, described below. In another embodiment, the coarse estimate may 

25 not be used, and instead, only the permutation parameter sets used to determined the 
refined values for the parameter variables. In yet another embodiment, a subset of the 
permutation parameter sets may be used to determined refined values for a subset of the 
N parameters. Thus, in various embodiments, at least a subset of the parameter sets 
comprising the coarse estimation and the permuted parameter sets may be used to 
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determine at least a subset of the N parameter variable values, i.e., the refined variable 
values. 

In the preferred embodiment, if there are N parameters estimated, then M 
permutations of the N parameters may be generated, where M is greater than or equal to 
5 N. As described below in more detail, M is greater than N so as to allow 
overdetermination of the N parameters. Each of the M permutations comprises a 
representation of a corresponding permutation waveform. In the preferred embodiment, 
where the pulse parameters are Op, t p , and (Op, N = 3. It has been determined that, due to 
diminishing returns of increased overdetermination, in the preferred embodiment, letting 

10 M = N+l provides for the greatest benefit in refining the N estimated parameters. 

In 208, in one embodiment, an inner or dot product, Qi, may be calculated 
between the received signal and the coarse estimate, as well as between the received 
signal and each of the M permutations of the estimated Gaussian Pulse parameters, i.e., 
between the received signal and each of the coarse estimate waveform and the M 

15 permuted Gaussian pulses or waveforms, generating M+l inner products. This 
calculation is facilitated by the fact that for Gaussian pulses there is a closed form, or 
analytic, solution for the inner product, as described below in the section titled Theory. 
Said another way, each of the plurality of parameter sets may be used in the equation for 
Q, providing a corresponding plurality of inner product equations. It should be noted 

20 that, as shown in the Theory section below, the closed form equation of the inner product 
includes terms for the estimation parameters, as well as variables for the as yet 
undetermined parameters of the received Gaussian pulse. In other embodiments, the 
coarse estimate may not be used to generate an inner product, and instead, the inner 
products may be generated from at least a subset of the permuted waveforms, or, 

25 alternately, the coarse estimate and a subset of the permuted waveforms may be used to 
generate the inner products. 

In 210, in one embodiment, the plurality of inner products, Qi, may be used to 
generate a corresponding plurality of linear equations. In other words, the linear 
equations may be derived from the inner product equations. Each of the linear equations 
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may be a function of a respective one of the estimation parameter sets and corresponding 
N parameter variables of the Gaussian pulse. Details of this process are described in 
greater detail below in the Theory section. In a preferred embodiment the plurality of 
inner products comprises M+l inner products, which may be used to generate a 
corresponding plurality of M+l linear equations. In other words, the M+l linear 
equations may be derived from the M+l inner product equations. Thus, where there are 
three Gaussian parameters (otp, t p , and tOp) to be determined for the received Gaussian 
pulse, there are three or more permutations for each estimate, giving four or more inner 
products, and therefore, four or more linear equations in the parameter variables. As 
mentioned above, in the preferred embodiment, when N = 3, then M = 3, and so there are 
M + 1 = 4 linear equations to solve the N = 3 variables. 

Finally, in 212, the plurality of linear equations derived in 210 may be solved to 
determine the parameters of the Gaussian pulse. It should be noted that in the case of 
three parameters or variables, a minimum of three linear equations may be solved to 
determine the variable values. However, by overdetermining the variable values through 
the use of more than the minimum number of linear equations, an "averaging" effect 
results, providing more accurate values for the parameters. In the preferred embodiment, 
the overdetermination may be accomplished by simply using one more linear equation 
than the number of variables. As mentioned above, the benefits of using even more linear 
equations for even greater overdetermination are marginal, providing diminishing returns 
for the number of calculations involved. Using the (over) determined parameters of the 
received Gaussian pulse, various analyses may be performed. For example, in the TDR 
system of Figure 1, a transmitted Gaussian pulse and one or more reflected pulses (from 
the DUT 110) may each be analyzed to characterize an interconnect discontinuity in the 
DUT 110. In one embodiment, the parameters may be used to determine the nature and 
magnitude of the discontinuity. For example, the characterized Gaussian pulses may be 
used to determine the transfer function of the discontinuity, or to calculate the impedance. 
As another example, the value of the time shift parameter t p for each pulse may be used to 
locate the interconnect discontinuity or flaw in the device or transmission medium. 
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In one embodiment, the CPU 101 of the TDR system of Figure 1 may execute the 
one or more software programs to detect and characterize the transmitted Gaussian pulse 
according to the method described above. In the embodiment where the signal also 
includes the one or more reflected pulses, the CPU 101 may execute the one or more 
5 software programs to detect and characterize the reflected pulses, as well, wherein the 
detection and characterization of the transmitted Gaussian pulse and the one or more 
reflected Gaussian pulses may be useable to perform TDR analysis on the DUT 110. 

Thus, by using various embodiments of the above-described method, a Gaussian 
10 pulse in a noisy signal may be detected and characterized. In one embodiment, the 
method may be used to distinguish between an original pulse and a plurality of reflected 
pulses in a Gaussian TDR system, and to characterize connection discontinuities in the 
DUT 110, as described above with reference to Figure 1. 

Figure 3 - Flowchart of a Method for Making an Estimation of a Gaussian Pulse 

Figure 3 is a flowchart of one embodiment of a method for making the estimation 
of parameters for the Gaussian pulse comprised in the received digital signal, as described 
above in 202 with reference to Figure 2. It should be noted that in various embodiments, 
various of the steps presented below may occur in a different order than shown, or may be 
20 omitted, as desired. In other embodiments, additional steps may be included, as well. In 
the preferred embodiment, the estimation may be made using the "zoom-in" approach, as 
is well known in the art, and which is described in detail below. At a high level, the 
zoom-in approach involves determining which of a predetermined set of Gaussian 
waveforms, also termed "Gaussian windows", most closely matches the Gaussian pulse in 
25 the received signal. Once the closest match is made, that Gaussian waveform's 
characteristic parameters may be used as the estimation of the Gaussian pulse. It should 
be noted that in the TDR system described above with reference to Figure 1, the method 
described may also be used to generate estimations of parameters for the one or more 
reflected Gaussian pulses. 
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The received signal may comprise a first sequence of values, as may be produced 
by the digitizer 106. As Figure 3 shows, in 302 a current area of interest of the received 
signal may be determined, i.e., a start position and a stop position in the first sequence of 
values comprised in the signal may be determined. The current area of interest may 
5 comprise a second sequence of values which includes at least a portion of the first 
sequence of values. 

In one embodiment, the initial current area of interest may comprise most or all of 
the received signal. In the preferred embodiment, the number of data points or values in 
the area of interest (second sequence of values) is a power of two, such as 256, 512, 1024, 
10 etc. The determined area of interest is the part of the received signal to be compared 
against a given Gaussian window. In the preferred embodiment, each determined area of 
interest (but the first) may be centered at the optimal position of the previous Gaussian 
window's peak. 

In 304, a current Gaussian window corresponding to the determined area of 
15 interest may be selected. In one embodiment, a Gaussian window may be a synthesized 
Gaussian waveform whose parameters are known, and which includes a number of data 
points which corresponds to, or is proportional to, the length of the determined area of 
interest described in 302. In other words, the Gaussian window may comprise a third 
sequence of values representing an elementary Gaussian function. In the preferred 
20 embodiment, the number of data points or values comprised in the selected Gaussian 
window is also a power of two, and is half that of the determined area of interest, such 
that if the initial area of interest has a length of 2048, then a Gaussian window of length 
1024 is selected. In other embodiments, the ratio of the lengths of the area of interest and 
the Gaussian window may be greater or less than two, 
25 In 306, a windowed FFT may be performed with the selected Gaussian window 

over the determined area of interest to generate a corresponding power spectrum. In one 
embodiment, the windowed FFT may be performed by logically aligning the Gaussian 
window at the start position of the area of interest such that each data point or value in the 
Gaussian window corresponds to a particular value in a portion of the area of interest, i.e., 
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wherein a sub-sequence of the values comprised in the area of interest aligns with the 
third sequence of values comprised in the Gaussian window. 

An element-wise product may be calculated between the Gaussian window (the 
third sequence of values) and the aligned portion of the area of interest (the sub-sequence 
5 of values), generating a product waveform. A Discrete Fourier Transform (DFT) may 
then be applied to the product waveform to determine a power spectrum of the product 
waveform. A peak frequency amplitude may be identified from the power spectrum. 
Once the power spectrum of the product waveform is determined, the Gaussian window 
may be shifted with respect to the area of interest, and the process repeated, generating 

10 another product waveform, and another power spectrum. In one embodiment, the 
Gaussian window may be shifted by one or more data points or values with respect to the 
area of interest data. The Gaussian window may then be shifted again, and so on, until 
the entire area of interest has been compared (via inner product and FFT) to the Gaussian 
window waveform. Thus, a peak frequency amplitude may be determined for each of the 

15 plurality of generated power spectra. The peak frequency amplitude for each generated 
power spectrum is a measure of the degree to which the Gaussian window/waveform at 
that position correlates with the waveform of the area of interest. It should be noted that 
the element-wise multiplication of the Gaussian window with the portion of the area of 
interest of the received signal, followed by the FFT of the resulting product is functionally 

20 equivalent to performing an inner product of the portion of the area of interest of the 
received signal with a modulated ideal Gaussian waveform (window). 

In 308, a maximum peak frequency amplitude may be determined from the 
plurality of power spectra generated in 306. In other words, the greatest peak value from 
all of the selected Gaussian window's power spectra may be determined, indicating 

25 maximal correlation of the waveform at the corresponding position. 

As Figure 3 shows, the process described in 302-308 may be repeated for each of 
the other Gaussian windows, generating corresponding maximum peak frequency 
amplitudes for each. In the preferred embodiment, at each iteration the length of both the 
area of interest and the Gaussian window is decreased by a factor of two. 
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Finally, in 310, the maximum peak frequency amplitudes determined in each 
iteration may be compared, the greatest maximum peak frequency amplitude determined 
from them, and the corresponding Gaussian window or waveform and corresponding 
position selected as the estimation of the Gaussian pulse. 

Thus, with each iteration the method "zooms in" by a factor of two, measuring the 
correlation between each Gaussian window or waveform at a multitude of relative 
positions with the area of interest waveform. The Gaussian waveform at the position 
with the highest correlation may be selected as the estimation of the Gaussian pulse, and 
its N parameters used in the permutation step described above in 206. 

As mentioned, in one embodiment, the method described above with reference to 
Figures 2 and 3 may be applied iteratively to a signal containing multiple Gaussian 
pulses, such as may be produced by the TDR system of Figure 1. Such as method is 
described below with reference to Figure 4. In this embodiment, each time a pulse is 
detected and characterized, the characterized pulse may be subtracted from the signal, 
leaving a residue containing any remaining pulses. This residue may then be used as the 
input signal, and the next pulse detected and characterized. Again, the characterized 
pulse may be subtracted from the signal, leaving another residual signal, and so on, until 
all pulses of interested have been detected and characterized. Useful analyses may then 
be conducted on the resulting characterized pulses to extract useful information regarding 
the DUT, such as the number, nature, and location of connection discontinuities therein. 

Figure 4 - A Method for Performing Time Domain Reflectometry Usin g Gaussian Pulses 
Figure 4 is a flowchart of one embodiment of a method for performing Time 
Domain Reflectometry (TDR) using Gaussian pulses. This particular embodiment relates 
to using TDR to locate a fault, such as in a transmission medium. The signal may 
comprise a sequence of values which includes a plurality of Gaussian pulses, as well as 
noise. It is assumed that a signal has been received, for example, from a memory 
medium, a digitizer, or any other source, and is in a form suitable for analysis. 
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As Figure 4 shows, in 402 a Gaussian pulse comprised in the signal may be 
characterized, such as by applying the methods described above with reference to Figure 
2 and Figure 3. In one embodiment, the characterization of the Gaussian pulse may 
include producing a set of parameter values which characterize the pulse, and from which 
5 an approximation of the Gaussian pulse may be generated. In one embodiment, the 
parameter values may be stored for later use. In another embodiment, the parameter 
values may be used to produce an approximation or construction of the Gaussian pulse, 
which itself may be stored. 

In 404, the Gaussian pulse may be subtracted from the signal. In one 
10 embodiment, subtracting the Gaussian pulse from the signal may include generating the 
approximation or reconstruction of the Gaussian pulse from the parameters of 402, then 
subtracting the approximated Gaussian pulse from the signal, i.e., numerically subtracting 
a sequence of values comprised in the approximated or reconstructed Gaussian pulse 
from the sequence of values comprised in the signal, resulting in a remainder signal, or 
15 residual. 

In 406, a determination may be made as to whether there are more pulses of 
interest in the residual. If there are more pulses in the residual, then steps 402 and 404 
are performed again, this time using the residual signal, and thereby producing a second 
residual. Again, a determination may be made as to whether there are more pulses of 

20 interest in the residual, and so on, repeating steps 402 and 404 until it is determined that 
there are no more pulses of interest in the final remainder or residual signal. 

Once it has been determined in 406 that there are no more pulses of interest in the 
residual, then in 408 the characterized pulses that correspond to an original pulse and one 
or more reflected pulses may be determined. In one embodiment, the plurality of 

25 parameter sets which correspond to each of the characterized Gaussian pulses may be 
sorted by one or more parameter values, such as by time shift, in order to clarify the 
nature and relationships of the characterized pulses. In one embodiment, the parameters 
may be used to calculate other parameters, such as amplitude, upon which the parameter 
sets may be sorted. The order of the sorted parameters may be used to ascertain which of 
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the pulses is the original pulse, and which are reflected pulses. For example, if the pulses 
are sorted by ascending time shift, then the top most pulse is the earliest received pulse, 
and thus may be assumed to be the original pulse. This conclusion may be confirmed by 
checking whether the top most pulse also has the greatest amplitude, as would be 

5 expected. Similar types of analyses may be performed on the remaining signals to 
determine a primary reflected pulse, and any other reflected pulses, as desired. 

Finally, in 410 the location of the fault in the medium may be determined. In the 
preferred embodiment, the location may be determined by comparing the time shift 
values of the original pulse and the primary, or first, reflected pulse, producing a time 

10 difference At. The time difference At combined with the known speed of the signal in the 
medium may be used to calculate the distance from the signal source to the fault, thereby 
locating the fault. It should be noted that the fault location is but one of a variety of 
features of the medium which may be investigated using the above method, such as the 
magnitude and/or nature of interconnect discontinuities, impedance, etc. It should also be 

15 noted that the TDR process described above is applicable not only to transmission media, 
but to any DUT. For example, the above method or similar methods may be used to 
analyze characteristics, e.g., analyze interconnect discontinuities, of signal traces or vias 
in a printed circuit board (PCB) or a printed wiring board (PWB). 

20 Figures 5A-5E - Example Signals with Gaussian Pulses 

Figures 5A-5E illustrate example signal traces which include one or more 
Gaussian pulses, including reflected pulses. In one embodiment, these signal traces may 
be produced by a system such as the TDR system described above with reference to 
Figure 1 , using the methods described above with reference to Figures 2-4. 

25 

Figure 5A - Gaussian Pulse and Reflection Pulses 

Figure 5A is an example of a signal which contains a Gaussian pulse and two 
reflection pulses. The largest pulse (left) 502A is the initial or original Gaussian pulse. 
The middle pulse and rightmost smallest pulse are reflected pulses 504B, where the 
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i. 



middle pulse is a primary reflected pulse, and the smallest is a secondary reflected pulse. 
As Figure 5A shows, the signal is extremely clean, i.e., signal noise is neglibible. This 
example is not representative of typical real world signals, where the presence of noise is 
almost inevitable. 

5 

Figure 5B - Gaussian Pulses with Refined Estimates and Residue 

Figure 5B is an example of a Gaussian pulse and a partially coincident reflected 
pulse. As Figure 5B shows, there are three traces overlayed on the chart; an original 
signal 506, which includes an initial Gaussian pulse 502B shown in the left three fourths 

10 of the trace, and a higher frequency reflected pulse 504B shown on the right. As 
described above, the method has been applied iteratively to the signal, generating a pulse 
estimation for each of the pulses, which together form an estimated signal 508. The 
estimated pulses (signal) 508 are subtracted from the original signal 506 to produce the 
low amplitude waveform, or residual, trace 510 as shown. The amplitude of the residual 

15 510 is a measure of the accuracy of the estimations, in that the greater the coincidence of 
the estimation pulses 508 with the original pulses 506, the smaller the difference, or 
residual 510. 

Figure 5C - Multiple Gaussian Pulses in Noisy Medium 

20 Figure 5C is an example of a signal 506B which includes many Gaussian pulses 

and which also includes substantial noise. As Figure 5C shows, a primary Gaussian pulse 
of amplitude ~. 4, centered at -.13 on the time axis is followed by a similar reflected pulse 
with approximately the same amplitude, located at -.185. Additionally, there are many 
pulses of unknown origin with amplitudes between 0 and .2, as may be seen. Figure 5D, 

25 described below, illustrates the result of applying the present method to such as noisy 
signal. 

Figure 5D - Refined Estimates of Multiple Gaussian Pulses in Noisy Medium 
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Figure 5D shows the results of iterative application of the present method to the 
signal 506B of Figure 5C, namely, the estimated or reconstructed signal 508B, Figure 5D 
is an exemplary illustration of the utility of the present method for extracting useful 
information from extremely noisy signals. As Figure 5D shows, the method has detected 
and characterized many (-100) Gaussian pulses in the signal. The estimated signal 508B 
may be subtracted from the original signal 506B of Figure 5C to produce a residue, as 
shown below with reference to Figure 5E. 

Figure 5E - Multiple Gaussian Pulses in Noisy Medium with Refined Estimates and 
Residue 

Figure 5E illustrates the original signal trace 506B from Figure 5C overlayed with 
the estimated or constructed signal trace 508B from Figure 5D. To clarify the matter, a 
portion of the signal trace has been magnified, illustrating two adjacent Gaussian pulses, 
512A and 512B, As may be seen, the difference between the original and estimated 
signals is also shown, in the form of the residual 510. As mentioned above, the accuracy 
of the estimate is reflected graphically in the residual, with a lower amplitude residual 
indicating greater accuracy. In this example, measuring the residual as a percentage of 
the original signal gives a value of -13%. 

As mentioned above, once the method has been applied to the signal, the 
quantified parameters generated for each detected pulse may be sorted by one or more of 
the parameters and analyzed. For example, the detected pulses may be sorted by 
amplitude, wherein the pulse with the highest amplitude may be designated as the 
primary, or original Gaussian pulse. Typically, the pulse with the second highest 
amplitude is the first reflected pulse. The parameters of these two pulses may then be 
compared to deduce characteristics of the reflection event or fault. For example, the time 
shift t p of the two pulses may be compared to determine the location of the fault, and so 
on. The information gleaned from these first two pulses may be useful in determining the 
relevance of each of the other pulses, allowing secondary reflection pulses to be 
distinguished from incidental pulses with the same or even greater amplitudes. 
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Theory 

The theoretical basis of the above-described method is based upon the adaptive 
Gabor expansion, which is a special case of the so-called Matching Pursuit (MP) as is 
well known in the art. It is assumed that the use of moving Gaussian windows to 
generate the estimate h P) o of the Gaussian pulse in the signal s p (t), is well understood, and 
so the estimate h p ,o with parameter set (a p?0 , the inverse of the Gaussian pulse variance, 
t p ,o, the Gaussian pulse time shift, and ci)p,o, the carrier or base frequency) is assumed as 
given. The following describes the method used to refine the estimate. 

The signal s p (t) may be considered to be comprised of two orthogonal 
components, or basis functions, which are the Gaussian pulse, and an orthogonal noise 
component, or residue. If the estimate h P) o is assumed to be close to an optimal or ideal 
Gaussian elementary function, then the following relation holds: 



(1) 



where v p is the component of the signal parallel to the optimal function, s p+ i is orthogonal 
to v p , and X denotes the correlation between h P) o and s p+ i. The term s p+ j is used to denote 
the orthogonal component because it is the residue left after removal of the parallel 
component, i.e., the characterized Gaussian pulse. This residue may then be examined 
for further pulses, and so is considered to be the "next" signal, s p+ ]. Because h P) o is close 
to v p , X may be considered negligible, thus the inner product between s p (t) and h P)0 (t) may 
be simplified to 



(2) 



= Q(ocp )0 , t P)0 > (%(>) 



(3) 
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Note that the inner product Q(Op,o, t p ,o, t%o) may be written in closed form for a 
Gaussian function: 



Q(Op <%m) = 



A. 



p p,m 



(4) 



While the set of parameters (Op, t p , (Op) is unknown and needed to estimate the 



th 

Gaussian pulse, the set of parameters (0Cp, m , t p>m , (0p ?m ) represents the m testing point. If 
it is assumed that 



Q(ccp, 0 , t p , 0 , %o) =\(s p (t) 9 h pi0 (t))\ = (s pJ h 



(5) 



then, in principle, the set of unknown parameters (Op, t p , (%) can be solved by a 
group of equations with testing points (<Xp fin , t p>m , 0)p, m ), for m = 0, 1,2,... 



Because both sides of equation (4) are positive, the following relation holds: 



(s h ) 

* P 9 <z Ptm t Ptin a Ptm f 



a a 

p p,m 

a +a 

P P>m 



2{a n -va A 

\ p p,m' 



p,m 



a +a 

p p.m 



(6) 



With elementary mathematical manipulation, equation (4) can be regrouped as 



ax + bco +ca -d 



(7) 
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where 



x - a at 

p,m p p 



and a, b, c, and d are known coefficients given by 



a — t — t 

p,m p,n* 



(8) 



b^o) -co n , 



p,m" p,m p,m 



p,n*P>n~p,n 



(9) 



(10) 



and 



d — a In 

p,m 



p,m p,tn P>tn 



p,n" p,tt^p,n 



(ID 



where <f derives from the fact that the inner product Q may be approximated by the inner 
product of the original signal and the elementary Gaussian function, h(t), as expressed in 
equation (2), above. 



By applying the set of following test points, 



(Op,0 

(Op,0 

(Op, 

(Op, 
(Op, 



t P ,o, %o) 

tp,l, <Dp,i) 
tp,l, <%-l) 
tp,-l» %l) 
tp,-h <0p r l) 



(12) 



(otp, tp, o)p) can be solved from the following linear system: 
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4 



A 



x 
a, 



a,. n - ln 



'jk,0 



(s h ) 2 



( S p ^K p0 t pl o> p _ 1 ) 



P ' "p.ffp.tPp.O 



) 



/ 2 2x 



where the matrix A is given by 



A = 



a 



p,0 



(s h ) 



a: 



p,0 



(^ 2 -^o 2 )" ln 



(5 Zi ) 



(s h 



> 



a 



^,0 - ^p.l 



(Vi"^oV ln 



2 W 



(5 A > 



P' S.o^Ao 



> 



a 



p,0 



P"l 



(f^ 2 -f Pi o 2 )-ln 



( S p ' K^^ea^ ) 



(13) 



(14) 



The set of parameters (ctp >0 , t p , 0 , c%o) denotes initial values which may be 
estimated by the "zoom-in" approach described above with reference to Figure 3. The 
10 remainder of the testing points may be determined by: 



t P ,n = t P ,o + nA t 



n = ±l 



(15) 



and 



%n = CDp,o + nA w n = ±1 



(16) 



15 
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The time increment A t is inversely proportional to the initial variance oCp t o. The 
larger the variance Op,o, the smaller the time increment A t , and vise versa. On the other 
hand, the frequency increment A® is proportional to the variance Op ; o. In other words, the 
larger the variance oc p? o, the larger the frequency increment A t0 , and vice versa. And so, by 
permuting the values of the coarse estimate parameters, multiple sets of variable values 
may be generated. These additional variable sets may each be used in the linear equation 
derived from the closed form inner product equation (4), thereby generating 
corresponding additional linear equations which may be used to solve for the refined 
variable set, e.g., (a p , t p , oOp). 

Although in principle, any optimization algorithm may be used to refine the 
estimate, such as the Newton-Raphson method, the approach described herein has been 
found to be far better in terms of speed and stability. The generic optimization algorithms 
often have convergence problems in multi-dimensional optimization applications. The 
success and efficiency of the method presented herein hinges largely on the validity of 
equation (5). The method fails when different Gaussian components of the signal are too 
close to be separated. 

As mentioned above, once the best match Gaussian window or elementary 
function has been determined, i.e., h p (t), the residual may be computed by 

s p+1 (t) = s p (t)-B p hp(t) 

where B p is the coefficient of the component corresponding to the selected best match 
Gaussian window. The method described herein may then be repeated to compute the 
next elementary function h p +i(t)[i]. The process may be repeated iteratively until ||s p+n (t)|| 
> e, where e is a predetermined error threshold. 
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Figures 6A-31B - Lab VIEW Virtual Instrument Set to Perform Gaussian Detection and 
Characterization 

Figures 6A-31B illustrate Virtual Instruments (Vis) which implement one 
5 embodiment of the present invention. The Vis comprise an example high-level 
implementation of one embodiment of the system and method described above. This 
embodiment was implemented in the National Instruments Lab VIEW development 
environment, and comprises a Lab VIEW Virtual Instrument (VI). Each VI has been 
implemented using the National Instruments Lab VIEW graphical programming language. 
10 Each VI is presented in two parts — a configuration panel (Figure nA), and a block 
diagram (Figure nB). 

Figures 6A and 6B - Adaptive Gabor Expansion VI 

Figures 6A and 6B illustrate a configuration panel and block diagram, 
15 respectively, for an Adaptive Gabor Expansion VI. The Adaptive Gabor Expansion VI 
uses weighted adaptive Gaussian functions to reconstruct the time waveform. The 
information of weights and Gaussian functions are computed by the Adaptive Gabor 
Transform VI. 



20 Figures 7A and 7B - Real Gaussian Pulse Generator VI 

Figures A and B illustrate a configuration panel and block diagram, respectively, 
for a Real Gaussian Pulse Generator VI The Real Gaussian Pulse Generator VI accepts 
characteristic parameters for a specified Gaussian pulse and generates the real-valued 
Gaussian pulse. 

25 

Figures 8 A and 8B - Ramp Pattern VI 
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Figures A and B illustrate a configuration panel and block diagram, respectively, 
for a Ramp Pattern VI. The Ramp Pattern VI generates an array containing a ramp 
pattern. 

Figures 9 A and 9B - Adaptive Gabor Transform VI 

Figures 9A and 9B illustrate a configuration panel and block diagram, 
respectively, for an Adaptive Gabor Transform VI. The Adaptive Gabor Transform VI 
performs an Adaptive Gabor Transform of an input sequence. 

Figures 10A and 10B - Multi-Scale STFT VI 

Figures 10A and 10B illustrate a configuration panel and block diagram, 
respectively, for a Multi-Scale STFT VL The Multi-Scale STFT VI performs a multi- 
scale windowed Fourier Transform of an input sequence. 

Figures 1 1 A and 1 IB - Unit Energy Gaussian VI 

Figures 11A and 11B illustrate a configuration panel and block diagram, 
respectively, for a Unit Energy Gaussian VL 

Figures 12A and 12B - Real FFT VI 

Figures 12 A and 12B illustrate a VI configuration panel and block diagram, 
respectively, for a Real FFT VI. The Real FFT VI computes a Real Fast Fourier 
Transform or a Real Discrete Fourier Transform of an input sequence. 

Figures 13A and 13B - Extension (real) VI 

Figures 13A and 13B illustrate a configuration panel and block diagram, 
respectively, for an Extension (real) VI. The Extension (real) VI generates extended data 
for an input sequence subject to specified constraints. 

Figures 14A and 14B ~ Linear Fit Coefficients VI 
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Figures 14A and 14B illustrate a configuration panel and block diagram, 
respectively, for a Linear Fit Coefficients VI The Linear Fit Coefficients VI finds the set 
of linear coefficients, i.e., slope and intercept, which describe the line that best represents 
an input data set. 

Figures 15 A and 15B - General Polynomial Fit VI 

Figures 15A and 15B illustrate a configuration panel and block diagram, 
respectively, for a General Polynomial Fit VI. The General Polynomial Fit VI finds the 
polynomial curve values and the set of Polynomial Fit Coefficients which describe the 
polynomial curve that best represents an input data set. 

Figures 16A and 16B - General LS Linear Fit VI 

Figures 16 A and 16B illustrate a configuration panel and block diagram, 
respectively, for a General LS Linear Fit VI The General LS Linear Fit VI finds the k- 
dimension linear curve values and the set of k-dimension linear fit coefficients which 
describe the k-dimension linear curve that best represents an input data set using the least- 
squares method. 

Figures 17 A and 17B - A x Vector VI 

Figures 17 A and 17B illustrate a configuration panel and block diagram, 
respectively, for an A x Vector VI. The A x Vector VI performs a multiplication of an 
input matrix and an input vector. 

Figures 18A and 18B - Mean Square Error VI 

Figures 18A and 18B illustrate a configuration panel and block diagram, 
respectively, for a Mean Square Error VI. The Mean Square Error VI computes the mean 
squared error of input sequences X and Y. 

Figures 19A and 19B - ID Linear Evaluation VI 
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Figures 19A and 19B illustrate a configuration panel and block diagram, 
respectively, for a ID Linear Evaluation VI The ID Linear Evaluation VI performs a 
linear evaluation of an input array. 

Figures 20 A and 20B - ID Polynomial Evaluation VI 

Figures 20A and 20B illustrate a configuration panel and block diagram, 
respectively, for a ID Polynomial Evaluation VI. The ID Polynomial Evaluation VI 
performs a polynomial evaluation of an input data set X using coefficients a. 

Figures 21A and 21B - Coarse Estimation VI 

Figures 21 A and 2 IB illustrate a configuration panel and block diagram, 
respectively, for a Coarse Estimation VI The Coarse Estimation VI generates a coarse 
estimation for an input sequence. The coarse estimation comprises a set of characteristic 
parameters which specify a Gaussian pulse. 

Figures 22 A and 22B - Refinement VI 

Figures 22A and 22B illustrate a configuration panel and block diagram, 
respectively, for a Refinement VI. The Refinement VI generates a refined Gaussian pulse 
estimation from a received coarse estimation. 

Figures 23A and 23B - Pseudolnverse VI 

Figures 23A and 23B illustrate a configuration panel and block diagram, 
respectively, for a Pseudolnverse VL The Pseudolnverse VI finds a pseudo-inverse 
matrix of a rectangular real matrix. 

Figures 24 A and 24B - Generate Matrix A for Gaussian Detection VI 

Figures 24A and 24B illustrate a configuration panel and block diagram, 
respectively, for a Generate Matrix A for Gaussian Detection VI. 
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Figures 25 A and 25B - Inner Product with Gaussian Pulse VI 

Figures 25A and 25B illustrate a configuration panel and block diagram, 
respectively, for an Inner Product with Gaussian Pulse VI. The Inner Product with 
Gaussian Pulse VI computes an inner product between an input sequence and a specified 
Gaussian pulse. 



Figures 26A and 26B - Complex Dot Product VI 

Figures 26A and 26B illustrate a configuration panel and block diagram, 
respectively, for a Complex Dot Product VI. The Complex Dot Product VI computes a 
10 complex dot product of complex vectors X and Y. 



Figures 27A and 27B - Complex Gaussian Pulse Generator VI 

Figures 27A and 27B illustrate a configuration panel and block diagram, 
respectively, for a Complex Gaussian Pulse Generator VI. The Complex Gaussian Pulse 
15 Generator VI generates a complex Gaussian pulse from specified characteristic 
parameters. 



Figures 28A and 28B - Generate Vector Y for Gaussian Detection VI 

Figures 27A and 27B illustrate a configuration panel and block diagram, 
20 respectively, for a Generate Vector Y for Gaussian Detection VI. 



Figures 29A and 29B - Updata STFT VI 

Figures 29A and 29B illustrate a configuration panel and block diagram, 
respectively, for an Updata STFT VI. 



Figures 30A and 30B - Dot Product VI 

Figures 30A and 30B illustrate a configuration panel and block diagram, 
respectively, for a Dot Product VI. The Dot Product VT computes the dot product of input 
vectors X and Y. 
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Figures 31 A and 3 IB - Co m pute Projection VI 

Figures 31A and 3 IB illustrate a configuration panel and block diagram, 
respectively, for a Compute Projection VI. The Compute Projection VI computes a 
projection of an input Gaussian pulse on a Gaussian Waveform basis vector. 

Figure 32- VI Hierarchy of the Adaptive Gabor Expansion VI 

Figure 32 illustrates the hierarchy of component Vis comprised in the Adaptive 
Gabor Expansion VI, described with reference to Figures 6A and 6B. As Figure 32 
shows, the Real Gaussian Pulse Generator VI described with reference to Figures 7A and 
7B, and the Ramp Pattern VI described with reference to Figures 8A and 8B are 
component VTs of the Adaptive Gabor Expansion VI. 

Figure 33- VI Hierarchy of the Adaptive Gabor Transform VI 

Figure 32 illustrates the hierarchy of component Vis comprised in the Adaptive 
Gabor Transform VI, described with reference to Figures 9 A and 9B. As Figure 33 
shows, the Adaptive Gabor Transform VI comprises a plurality of component Vis, 
specifically, the Vis described above with reference to Figures 10A-31B. It should be 
noted that some of the component Vis are sub-components, in that they are component 
Vis of a VI which is itself a component VI of the Adaptive Gabor Transform VI, and so 
on. In other words, the VI hierarchy has a depth greater than two. 

Example Applications 

In addition to the above-mentioned uses in the field of Time Domain 
Reflectometry, various embodiments of the above described system and method may be 
useful in analyzing an electroencephalogram (EEG), which is a record of the electrical 
activity of the brain made from electrodes placed on the scalp. The signal is extremely 
complex and reveals both rhythmical and transient features. One such feature is referred 
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to as a "sleep spindle". The definition given in a classical reference states the following: 
"The presence of sleep spindle should not be defined unless it is at least 0.5s duration, 
i.e., one should be able to count six or seven distinct waves within the half-second period. 
(...) The term should be used only to describe activity between 12 and 14 cps (cycles per 
5 second)." 

Although intended for standardization of visual detection, the above definition is 
almost directly translatable into the time-frequency parameters of adaptive Gaussian 
functions. The shape of sleep spindles is so similar to Gaussian functions that one can 
assume a one-to-one relationship between sleep spindles and adaptive Gaussian functions 

10 fitted to EEG. Thus, embodiments of the present method may be used to detect and 
characterize sleep spindles in an EEG signal. Compared to other time-frequency 
representation schemes, embodiments of the present invention offer extra advantages: 
high resolution, fast computation, local adaptability to transient structures, and 
compatibility of its time -frequency parameters with definitions of EEG structures used in 

15 the traditional, visual analysis. The detection of sleep spindles is but one example of the 
use of the present method to analyze EEG signals. 

Dr. Durka of the Laboratory of Medical Physics, Warsaw University, has said, 
"The computer revolution in EEG analysis still consists mainly of the fact that what was 
on paper now is watched at the monitor screen, mainly because the new signal processing 

20 methods were incompatible with knowledge about visual EEG analysis. Matching 
Pursuit (more specifically, the adaptive Gabor expansion) corresponds directly to the 
traditional way of analyzing EEG, namely visual analysis. For over 40 years of clinical 
practice, structures present in EEG being described in terms like 'cycles per second' or 
'time width, 1 are directly interpretable in terms of time-frequency parameters. Matching 

25 Pursuit is one of the first methods that allow for direct utilization of huge amounts of 
knowledge collected over years of clinical EEG." 

Thus, the above-described method may be used to detect, characterize, and 
analyze EEG data in a manner which is substantially more efficient and quantitatively 
more accurate than prior methods based on subjective visual inspection. It should be 
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noted that the above applications are merely illustrative, and are not meant to limit the 
domains of use of various embodiments of the present invention. 

Other applications of the described systems and methods include financial, stocks 
and bonds, and currency analysis, as well as medical and scientific data characterization 
5 and analysis. In other words, various embodiments of the present invention may be used 
to analyze any data sets or signals which include Gaussian waveforms, or approximations 
thereof. 

Thus, the systems and methods presented above provide a fast and efficient 
approach to detecting and characterizing modulated Gaussian pulses in a noisy signal. 
10 Additionally, the systems and methods described provide an improved approach for 
performing Time Domain Reflectometry to analyze a Device Under Test using Gaussian 
pulses. Finally, the use of Gaussian pulses allows the above systems and methods to be 
implemented with a personal computer and relatively inexpensive components, granting 
substantial savings over prior systems (as much as 90%), 

15 

While the present invention has been described with reference to particular 
embodiments, it will be understood that the embodiments are illustrative and that the 
invention scope is not so limited. Any variations, modifications, additions, and 
improvements to the embodiments described are possible. These variations, 
20 modifications, additions, and improvements may fall within the scope of the inventions as 
detailed within the following claims. 
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